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This paper presented a compressive sensing (CS) calibration technique that 
simultaneously measures all the element excitations of a 2x2 active phased 
antenna array (APAA) instantaneously. The power patterns of the 
rectangular 2x2 APAA element are synchronized to orthogonally cross 
multiplicative sub-array system applying compressive sensing in achieving 
an array thinning along two 1-D sub-arrays for a fixed steered beam 
radiation. The far field radiation pattern direction of the 2x2 APAA is 
examined considering the mutual effect, complex excitation of the amplitude 
and phase of the antenna elements placed on the x and y axes. The unwanted 
array excitation with errors are calibrated using minimization vector as a 
standard basic pursuit problem in compressive sensing technique for the 2x2 
APAA calibration. The elements amplitude and phase shift variation are 
compared with reference state 0 in order to evaluate the faulty error array 
elements. The performance evaluation of this proposed technique is 


measured and demonstrated to validate the effectiveness of the proposed 
antenna calibration technique. 
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1. INTRODUCTION 

Now adays, identification of failures in antenna array systems is theoretically and practically studied 
for efficient signal propagation in radar communications, satellite, massive multiple-input and multiple-output 
(MIMO), wireless communications and military applications [1]. Various advanced technologies have been 
developed and investigated over the past few years that have the potential to increase the capacity of a 
wireless communications system. In order to solve the problem of defficent elements in an antenna array 
system, many approaches have been proposed in literatures including; back propagation technique, matrix 
technique, neural networks, case based reasoning, genetic algorithms, source reconstruction technique and 
signal processing based algorithm [2]-[9]. 

Compressive sensing has been applied in various approaches to failures in some elements in an arrays 
under test (AUTs) to detect the faulty elements in a phased or active array system [10]-[14]. To resolve such a 
problem, initial knowledge of the failure-free antenna array radiation pattern, amplitude, and phase is required 
and it is referred to as reference state or state 0, in examining the faulty array elements [15]. Hence, a fast 
testing, diagnosis, and calibration of antenna array systems are essential. 

Calibration starts with adequate analysis of uncertainties to define the error tolerance, selection of 
the appropriate calibration, and measurement technique. Considering the AUT, error factors can be classified 
as dynamic errors with time difference, amplitude, phase shift variation and static errors which includes the 
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manufacturing tolerances of the antenna array system [16]. Considering the categorization of the calibration 
process, an on-site calibration procedure is mostly implemented to compensate the error influence of the site 
or to update the coefficients because of the frequency change for transmission. 

Calibration is important in deriving an adequate antenna array propagation performance by 
compensating the errors in the system. Generally, there are two classes of calibration techniques namely; 
single element and multiple elements measurement techniques [17]. Therefore, different approaches have 
been applied in literature which includes; the phase toggling technique [14], rotating element electric field 
vector (REV) technique [18] and four orthogonal phase states technique [19]. A recursive matrix-forming 
approach for the matrix selection of a multi-element phased array calibration [17], and a modified matrix 
technique for phased array calibration by resolving the linear equations are investigated in this paper [20], 
phase-encoding matrix modified from the finite phase-encoding setting accuracy for active phased-array 
antennas for synthetic aperture radar (SAR) [21], fast Fourier transform (FFT) and harmonic characteristic 
analysis for parallel calibration of phased array [22], and Fourier coefficients of one array-power response 
measured with the conventional rotating element electric field vector (REV) technique [23]. 

Recently, compressive sensing (CS) theory is evolving in solving linear equations relating to array 
diagnosis for evaluating the amplitude and phase of a near and far field measurements [11]-[13]. 
Consequently, a compressive sensing approach is being applied to detect faulty radiators in planar antenna 
arrays using a reduced set of phaseless far-field measurements, demonstrating that the nonlinear problem can 
be solved as long as the percentage of faulty elements is low and a smart measurement strategy is used. 
Numerical simulations and experimental results, have confirmed the good performance of the CS technique 
in terms of probability of fault detection [24]. 

In this paper, 1; compressive sensing technique is applied to diagnosis and calibrate a faulty error 
elements in 2D active phased array antenna (APAA) system. The faulty antenna element errors in the APAA 
casues a slight amplitude and phase shift variation which will be calibrated with a reference state 0 to 
evaluate the faulty elements. The unwanted array excitation with errors are calibrated using minimization 
vector (w) as a standard basic pursuit problem in compressive sensing approach for the 2x2 APAA 
calibration. In the following section, the research theory is presented in section 2, the methodology is presented 
in section 3, results and discussions is presented in section 4 and the conclusion follows in section 5. 


2. THEORY 

A 2D APAA with dimension P X Q AUT elements is uniformly distributed in a square equispaced 
grid of the p, elements in the x axis and qy elements on the y axis. To avoid grating lobes and high side 
lobes, r, and r, should be <//2. The array factor (æ) of the isotropic APAA is described in the (1) as: 


2T 


A TS a 
x (Ø ,) = LEZ Any eile x sin8cosh+Dpx ) ELA 


1, 2TrY 
e1 z sinĝðcosh +Dpy) (1) 


qx 
Apx = antenna element amplitude excitation along the x axis 
Agy = antenna element amplitude excitation along the y axis 
Dpx = antenna element phase shift along the x axis 

Dgy = antenna element phase shift along the y axis 

T = spacing between px elements along the x axis 

Ty = spacing between q, elements along the y axis 
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The direction of the APAA far field radiation pattern is (0,@), which relates to amplitude variation, phase 
shift variation, and radiation pattern. In a typical measurement setup, Figure 1 shows an APAA calibration 
system with each antenna element connected to a digital phase shifter controller. Because of the phase 
shifter’s amplitude and phase inaccuracies, higher order Fourier series are usually used to quantify array 
power responses. Consequently, this research focus on an adequate and efficient algorithm applying 
reconstruction of 2D signals analysis shown in the (1) to (3) on the missing samples. The 2D FFT of an 
active phased array antenna is expressed in the (3). Where, h = wave number. 


h== (4) 
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Figure 1. An active phased array antenna measurement calibration schematic 


The mutual effect, complex excitation of the amplitude, and phase shift variation of the antenna elements 
with respect to the x and y axes are evaluated when determining the far field radiation pattern direction of the 
P x Q APAA. The P, and Q, are sample along the x and y directions respectively for the 2x2 APAA as 
shown in the (7). The signal analysis is examined on a full set of G where: 


ni = hy — hyi; m; = hy — hyi (5) 


The hy; and h}; represents frequency of the ith component in the 2x2 APAA. 


G = {g(x,y) : xe{1,.., PL Y{L > Qy}} (6) 
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(7) 


Consequently, the error emmanting from the amplitude, phase shift and mutual coupling variations 
affects the radiation pattern of the 2x2 APAA. These errors are considered as the missing samples in the 2D 
spectral domain analysis. The zero-value expressed in the (8) act as an addition of noise to the signal, and the 
error is given as: 


_ (~g (x%y),f or (x,y) missing 
E(x, y) = lo,for (x,y) available (8) 


The noise is assumed to be zero-valued at the position of available signal samples, whereas the noise 
neutralizes the signal samples setting to zero at the position of missing samples. The available signal samples 
are expressed in the 2D FFT as: 


S(hy, Ry) = DD, (GG) + ey} (9) 


To analyze the error variations in the signal, the following cases are considered in the (10) to (13). 
- Ifh, = hzi, hy = hyi; the assumption is expressed as: 
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- Ifh, # hyi hy + hyi; the assumption is expressed as: 


yi» 
E{S(h,, hy)} = 0, 
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The noise affecting the P x Q APAA is expressed in s(h,, hy) as: 
o? {S(hy # hyi, hy # hyi)} = Pr Qype E{G (x, y)Q* (x, y)} 
* -1 
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Signal components (useful signals) and noise components can be distinguished using estimated 
noise parameters such as amplitude, phase shift, and mutual coupling fluctuations. As a result, the level of 
noise can be estimated, and the position of actual signal components above the noise can be detected, and the 
impacted elements can then be calibrated to the original signal. The entire variance of noise in the 2D FFT 
due to amplitude, phase shifter errors, or external errors (noise) on the elements can be stated in the (13) as: 


— çh 2 PxQy— PxkQyk 
g? = i=1 Zik Paye poa (13) 


3. RESEARCH METHODOLOGY 

In the previous section, we demonstrared how to analyze the noise variance in the 2D APAA caused 
by the amplitude, phase shift, mutual coupling variations and external noise disturbances. To attain a reliable 
array diagnosis technique that detects the defective elements position with minimal diagnosis time while 
calibrating the faulty elements to the original signal known as the reference state or state O at the receiver, 
as shown in Figure 2, it is necessary to investigate the effect of noise disturbance in the APAA. Calibration 
technique estimate and compensate antenna array errors, as well as associated errors and their location. 


Beam Steering 
Command (9, ) 


Amplitude Cmsreniè Sonah Phase 
Calculation ae me Calculation 
and Weight Vector 
(A) lo) 


Figure 2. Compressive sensing and vector weight technique for an active phased array antenna calibration 


3.1. Fault detection and calibration at the receiver 
3.1.1. Description of the method 

The vector of measurements y containing incorrect Pyg Qyg elments out of the original P,Q, signal is 
expressed as: 


f =s(Q) (14) 


The analysis of the Fourier transform of the signal is represented in the 2D FFT matrix form as 
D (Pr Qyx © P.Qy), this is extracted from the KRONECK PRODUCT of two FFT matrices as expressed in 
the (15). 


D = FFT PrQyx Q FFT P,Q, (15) 


The unobstructed APAA output signal f(t) is derived by multiplying the complex weight vector (u) and 
signal vector D(t) as shown in Figure 2 to give. 
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The complex factors Z„ in each branch models the amplitude, phase shift, and mutual coupling effects or 
external sounds deviation of the signal. The deviations in complex notation are described by the diagonal 
entries Z,,,n €{1, ... N}. Hence, the disturbed output signal is expressed as: 


FŒ = UT. DŒ) = U ... Uy. (17) 


Z,D,(t) 


f = uT = U'ZD(t) (18) 


ZyDy(t) 
The disturbed weight vector U7 is given as: 
F(© = UTZD(t) = (UZ, ... UyZyn) PD HD) (19) 


The angular frequency (w) and transmission angle (0) controls the P x Q APAA calibration by the steering 
vector of the array biased on the weight vector as shown in the (20). 


: 0 r 0 
We oe A E A =, WZ,COS— WZN-1cos— 
Ucar = Ucar w, 0) = }1, e? c, ye Cc (20) 


Where, c = speed of propagation. 
Applying the calibrated weight vector compensates the deviation in the array branches in the beamforming as: 


f(t) = OF, (w, 0) ZD(t) (21) 


3.2. APAA error detection and calibration technique 

Figure 3(a) shows the components of the APAA which includes, the antenna array, power divider 
and the 6 bits phase shifter. The experimental procedure of this research study on the the 2x2 APAA is 
carried out in Figure 3(b) with the antenna elements radiating in free space and the amplitude and phase shifts 
of this AUT is measured at far field without failure and expressed as the reference state or state 0. The 2x2 
antenna reflection coefficetion is investigated and measured using a vector network analyzer (VNA), 
showing the operating frequency at 3.5 GHz as displayed in Figure 3(c). Hence, modifying the array factor 
(a) in +z hemisphere is given as: 


A = sin 0 cos Q 
B = sin 0 sing 


x (A,B) = D (US, .Ugy) (22) 


Consdering the complex weight U distributed along the x and y directions of the P x Q APAA, the power 
pattern is expressed as: 


P =x (A,B) x « (A, B)* (23) 
(x) = element by element multiplication 


(*) = complex conjugate operation 
Therefore, applying the distributive property on the (23), it can be expressed: 


P=(D (Ups) x (D (Upx))") .(D (Ugy) x (D Uay))") (24) 


On the (24) illustrates that a power pattern is produced by substituting multiplication with an auto-convolution 
(®) of the complex weighting elements. 
The rectangular P Xx Q APAA power pattern employes two orthogonal 2D arrays written as: 


P = (DUR ® Uf))*).((PUoy 8 Uoy))) 25) 
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The combining of the two 1-D element patterns necessitates the use of a separate post processing 
unit at the receiver in order to obtain the necessary spatial power pattern function. The synthesizing pattern 
evaluates the complex weights (Umx, Uny) and distance between elements (pdx, qdy) along x and y axes. 
In the compressive sensing technique for the APAA calibration system, minimization of vector U is a typical 
basic pursuit obstacle. The expression below in the (26) and (27) represents the 2-D power pattern along the 0 
and ¢ plane. 


[POT cme = Ae x ilie (26) 
[P(®)ixey = [UylHaoy * [Exloy (27) 


Where: 

P(@) = power pattern splits along 0 Plane 

P(@) = power pattern splits along @ Plane 

Ux = x plane complex excitation vector in the P x Q APAA 
Uy = y plane complex excitation vector in the P x Q APAA 
L = matrix encoding the steering vectors 


4. RESULTS AND DISCUSSION 

This section presents a practical example for calibrating a 2x2 APAA as shown in Figure 3(a). 
Before imposing amplitude, phase shift, and mutual coupling changes or external noises deviations on the 
array, the four planar elements were simulated uniformly to determine their reference state 0. The amount of 
computation for the proposed method is examined and the factors influencing the calibration accuracy are 
analyzed qualitatively. 


4.1. 2x2 APAA design and fabrication 

A 2x2 APAA has been designed and fabricated using a 0.018 mm thick Rogers RO4350 Substrates. 
To achieve a sufficient coupling condition for the rectangular geometry of the array, as illustrated in Figure 3(b), 
the feeding ports are carefully determined along the four ports power divider used to achieve isolation between 
the input ports while maintaining a matched condition of all ports. Figure 3(c) illustrates the operating 
frequency of each of the four elements in the APAA with a maximum gain of 12.50 dB at reference state 0. 

Therefore, each element in the array is coupled to a 6-bits phase shifter module for a microwave 
network module to offer a configurable phase shift of the radio frequency (RF) signal. Figure 3(b) shows how 
the 6-bit phase shifter employed in this work is used in phased arrays systems for high-frequency RF 
applications. These phase shifters are coupled to a phase shifter controller, which has a four-way phase shift 
output and is regulated by a four-way dual in-line package (DIP) switch to measure various states while in 
the high and low states. Experimental studies using the Rohde and Schwarz network analyzer as shown in 
Figure 3(c) illustrates the return loss measurements of the APAA at 3.5 GHz. 


(a) (b) (c) 


Figure 3. Photos of the experiment system: (a) 2x2 array, power divider and 6-bits phase shifter; (b) 2x2 planar 
array experiment system in anechoic chamber; and (c) 2x2 array reflection coefficient measurement at 3.5 GHz 


Many array far-fields are measured, illustrating the excitation error as the array is uniformly excited 
considering the error variation. The 2 FFT is utilized to analyze the measured signals of the 2x2 APAA 
system, as shown in section 2. As the antenna elements are energized, these data signals are saved. The signal 
errors are detected using the scattering effect of the power divider in the array system, as well as faults caused 
by amplitude and phase shifter errors were evaluated using the 1, compressive sensing method, as shown in 
section 3. As a result, if the variation between the most recent calculated excitation error and the best among 


A compressed sensing approach for antenna array calibration (Akinwale Oluwaseyi Fadamiro) 


898 o ISSN: 1693-6930 


the previously calculated excitation errors is less than a target level, such as 10%, the average convergence 
error is considered the array excitation error under uniform array excitation, and the calibration procedure 
ends. The comparison of the genuine excitations and the calculated l} compressive sensing yields the 
calibration errors of this experimental calibration technique. Although, the calibration time increases due to 
the complexity of system but the accuracy and effectiveness of the method compensates for time estimate of 
the element under calibration (EUC) [25]. Figure 4(a) compares the pre-calibration, post-calibration, and 
normalized radiation patterns with no excitation errors. In the range of the sidelobes and nulls of the radiation 
patterns, the post-calibration pattern is substantially closer to the ideal pattern than the pre-calibration pattern. 
Because the sidelobes and nulls are particularly sensitive to the array excitation, the agreement between the 
post-calibration and the ideal radiation pattern suggests that the post-calibration excitations are very close to 
the ideal excitations. Figure 4(b) and Figure 4(c) compares the pre-calibration, post-calibration, and 
normalized radiation patterns with no excitation errors of the amplitude and phase excitation respectively. 
The undesirable array excitation with inaccuracies must be calibrated since the actual array excitation affects 
calibration accuracy. In practice, we can diagnose and calibrate antenna arrays under various array 
excitations and select the most appropriate technique. Finally, the proposed calibration technique is validated 
experimentally. 


Normalized Radiation Pattern (dB) 
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Figure 4. The idea, pre-calibration and post-calibration comparative analysis the 2x2 APAA measurement: 
(a) radiation pattern, (b) amplitude, and (c) phase 
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5. CONCLUSION 

This paper proposes a technique that can calibrate a 2x2 APAA simultaneously. The signal 
measurements of the four elements in the APAA are analyzed with the reference state 0 where no errors exist 
in the system. The key factor of this research study is to identify the problem from amplitude, phase shift, and 
mutual coupling fluctuations or external noise deviation errors on the APAA and use comprehensive sensing 
to model to resolve it. The toggling of the 6-bit phase shifters attached to each antenna element during the 
measurement operation determines the accuracy, complexity, and hardware requirements of this approach. 
To accomplish acceptable calibration, the phase settings are calibrated using 2D FFT signal data processing 
and basic pursuit compressive sensing combined with the vector weight (u). The efficiency of the suggested 
APAA calibration technique is demonstrated in an experimental evaluation of a far field 2x2 active phased 
array microstrip antenna array. It’s worth noting that these findings apply to radar communications, satellite 
communications, massive MIMO, and military applications. 
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